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Abstract 

An identity by Chaundy and Bullard writes 1/(1 — x) n (n — 1,2,...) as a sum of two 
truncated binomial series. This identity was rediscovered many times. Notably, a special 
case was rediscovered by I. Daubechies, while she was setting up the theory of wavelets of 
compact support. We discuss or survey many different proofs of the identity, and also its 
relationship with Gaufi hypergeometric series. We also consider the extension to complex 
values of the two parameters which occur as summation bounds. The paper concludes with 
a discussion of a multivariable analogue of the identity, which was first given by Damjanovic, 
Klamkin and Ruehr. We give the relationship with Lauricella hypergeometric functions and 
corresponding PDE's. The paper ends with a new proof of the multivariable case by splitting 
up Dirichlct's multivariable beta integral. 



1 Introduction 

Chaundy and Bullard noted "in passing" the identity 



i = (i - *r 1 f; ( n + ^ x k + *)*' 

k — k — 

as a side result in their 1960 paper John Smith's problem, see [31 p. 256]. Here m,n are nonneg- 
ative integers. Formula (|1.1|) can be written more succinctly as 

Pm,n( X ) + Pn,m(l ~ x) = 1, (1.2) 

where 

Pm , n ( x ) .= (i - zr 1 jr( n+ k k )x k = (i- f; (1.3) 

k=0 ^ ' k=0 

ja(a + l)...{a + k-l) if k = 1, 2, . . . , 

{a)k:= i if* = o, (L4) 



and 



is the Pochhammer symbol. 

The Chaundy-Bullard identity (jl.ip was rediscovered (partially or completely) many times: 



• In 1971 Herrmann [9J interpreted p m ,n(x) (see (|1.3p ) as the polynomial of degree m + n + 1 
which has a zero of order n + 1 at x = 1 and such that 1 — p mn (x) has a zero of order 
m + 1 at x = 0. He proved this by induction with respect to n (although we think that he 
meant induction with respect to m + n). Essentially, although not explicitly given in [9], 
Herrmann's result implies the identity (jl,2|) . 

• The identity (|1.1|) was proposed in 1975 for the Canadian Mathematical Olympiad (but 
not used there). Next it was proposed in 1976 for the problem section of Crux Mathe- 
maticorum. A proof by induction by Kleiman was given there [15] in 1977. The same 
identity was also proposed in 1977 for the elementary problem section in the American 
Mathematical Monthly by Burman. The Monthly [16J gave two solutions in 1979, one 
probabilistic proof by Schmitt and one using partial fractions by Jagers. Much later, in 
1992 in the Monthly [18] the probabilistic proof was implicit in the solution of a problem 
about the longest expected world series, posed in 1990 by Schuster. In 1997 the Monthly 
|19] had a follow-up with some non-probabilistic proofs. 

• A two- variable analogue of the identity (jl.lh was proposed in 1985 for the problem section 
in SIAM Review by Klamkin & Ruehr. In 1986 Bosch & Steutel gave in this journal [T7] 
a probabilistic proof as solution. In [17] it was also observed by Damjanovic, Klamkin & 
Ruehr that there is an n- variable generalization of the identity: 

5> E - E s ki , ai (/ y," ' ' + k * nh ^ . . . xt = i (x 1 + ---+x n = i). (i.5) 

i=l fc 1= o fc n =0 1 n ' 

They gave a proof by generating functions. A probabilistic proof was also indicated. 

• In 1988 Daubechies [U Lemma 4.4], see also [5j (6.1.7), (6.1.12)], rediscovered the case 
m = n of (jl.ip . This identity was a crucial step for her in order to arrive at the form of 
the function mo(£) which is associated with the wavelets of compact support named after 
her. Her proof in [5] was essentially the same as Jagers' proof in [16] , but she referred to 
Bezout's identitity. Next Zeilberger [22] in 1993 gave a probabilistic proof of Daubechies' 
case m = n of (jl.ip and, unaware of [17], he stated the case a\ = . . . = a n of (jl.5p and 
indicated a probabilistic proof. 

• Multiplication of both sides of (jl.ip by (1 — x)~ n_1 gives 

(i _ x yn-i = jr + x m+i (m ?, 1)fc (i - x) k ~ n -\ (i.6) 

fc=0 ' k=0 

This identity is the case m = of the identity at the end of section 8 in Vidunas [21] , where 
it is given as a three-term identity for three Gaufi hypergeometric functions satisfying the 
same hypergeometric differential equation in the most degenerate case (trivial monodromy 
group) . 

• It was also essentially identity (jl.6p which was rediscovered by Pieter de Jong (Nether- 
lands), who is studying the mathematical foundations of architecture. It was by his com- 
munication to the first author in 2007 that we first became aware of this identity. See also 
de Jong's manuscript [TP] . 
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It is not without precedents in mathematics, in particular in special function theory, that a 
relatively elementary result is rediscovered and published many times. We think that for (jl.ip 
the elegance and unexpectedness of the identity arose people's interest again and again, as it did 
with us. Why then spend again a publication on it? First it seems useful to survey all earlier 
(as far as we know now) occurrences and approaches. Second, we can offer some approaches 
which did not yet occur, notably the approach by splitting up the beta integral and the context 
of the Gaufi hyper geometric function. Third the possible generalizations are interesting. As 
we mentioned, the n-variable generalizations already occurred, but we can offer yet unexplored 
aspects of it. There are also various analogues and generalizations of (jl.ip in the g-case, which 
we will present in a forthcoming paper. 

The following sections present or survey many different proofs of (jl.ip . The proof in section [2] 
is the original proof by Chaundy and Bullard [3J, and its slight variations by Daubechies §6.1] 
and Jagers [16] are also discussed there. The proof in section [3] is by induction with respect to 
m + n. The proof in section H] is by repeated differentiation of (j4.1j) (a method suggested by 
Pieter de Jong in an earlier version of [TO]). The proof in section [5] uses generating functions. 
It was communicated to us by Helmut Prodinger and it is also a one-variable specialization of 
a proof in [T7]. Another proof of combinatorial flavour in section [6] uses weighted lattice paths, 
and becomes by specialization a probabilistic proof (different in formulation but in essence the 
same as many earlier proofs which appeared). Section [7] gives a proof by splitting up the beta 
integral. In section[8]we consider the extension of (jl.ip to complex values of m, n. In section[9]we 
observe that the three terms in the identity (|1.6p all solve a very special (degenerate) case of the 
hypergeometric differential equation. Next we obtain (jl.6p as a limit case of a more general three- 
term identity for hypergeometric functions, where the three terms all solve a hypergeometric 
differential equation. 

Section [TO] starts our discussion of the multivariable analogue (jl.5p . which was first obtained 
in [T7]. A connection with Appell and Lauricella hypergeometric functions is made. In section 
[TT1 a partial differential equation satisfied by all terms in this multivariable analogue is given. 
Finally we give in section[12]a proof of (jl.5p by splitting up Dirichlet's multivariable beta integral 
(generalizing the approach in section [7]) . 

Acknowledgements We thank Pieter de Jong for showing us his identity and suggesting the 
proof in section 01 Peter Paule for suggesting the approach with a multivariable beta integral in 
section [T2l and Mike Keane, Oleg Ogievetsky and Helmut Prodinger for helpful remarks. 

The second author was partly supported by FWF Austrian Science Fund grants P17563-N13, 
and S9607 (the latter is part of the Austrian National Research Network "Analytic Combina- 
torics and Probabilistic Number Theory"). 

Notation The Gaufi hypergeometric series (see [1] Ch. 2]) is defined by 

2 F 1 h b ;z) : =f;M^|i^ (z,a,b,ceC, \z\ < 1, c £ {0, -1, -2, . . .}), (1.7) 

V c / fc=Q \ c )k 

where the Pochhammer symbol is given by (jl.4p . In the terminating case we have 
2 F x (~ n ^-z\ :=E ( ~" )fc ffi* * fc M'CGC, n = 0,l,2,..., c + 0, -1, . . . , -n + 1). (1.8) 
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2 Chaundy &; Bullard's original proof 

Fix m and n. By the binomial theorem we have 

(x + y) m+n+1 = y n+1 P m , n (x, y) + x m+1 P n , m (y, x), (2.1) 

where 

P m ^y):=f^( m + n + l )x k y-K (2.2) 
fc=o ^ ' 
is a homogeneous polynomial of degree m. Put y := 1 — x. Then 

1 = (1 _ x ) n+l P mtn (x, l-x)+ x m+1 P nim (l -x,x), (2.3) 

and multiplication by (1 — x)~ n ~ l yields 

(1 - x)- n - 1 = P m , n {x, l-x) + x m+l (l - xy n - l P n)Tn {l -x,x). (2.4) 

Expand both sides of (|2.4p as a power series in x, convergent for \x\ < 1. Then P m ^ n (x, 1 — x) is a 
polynomial of degree < m in x and all terms in the power series of x m+1 (l — x)~ n ~ 1 P n<m (l — x, x) 
have degree > m + 1. Hence P mjn (x, 1 — x) equals the power series of (1 — x)~ n ~ l truncated 
after the term of x m , i.e., 

P m , n (x,l-x) = ^^±^x fe . (2.5) 

fc=0 

Then substitution of (|2.5p in (|2.3p proves (jl.ip . and its homogeneous form 

(x + y)^^ 1 = y^ 1 Y {n + 1)k x k (x + y) m ~ k + x m+1 V ^±^y k (x + y)"" fc . (2.6) 
k=0 k=0 

Note that, conversely, implies (|2.5|) . i.e., the equality 

£ ( m + k + x ) xfc(i - x)m_fc = e ( n t fc ) (2 - 7) 

fc=0 ^ ' k=0 ^ ' 

Indeed, compare (jl.6p with (|2.4p . It was essentially this identity (|2.7p which was also stated by 
Guenther [14J in a comment to the solution of a problem in the Monthly. He pointed out many 
relationships of this identity with the binomial and negative binomial distribution, including a 
probabilistic proof. 



Remark 2.1. We can rewrite (I2.7D as 

( n + 1 )kh_ n ^» {-m-n- l) k ( x 



fc! ' ^ k\ \x-l 

k=0 k=0 



In terms of terminating Gaufi hypergeometric series (jl.8p this can be written as 

/— m,n+l \ . „ f—m,—m — n — l x \ 

2Fi[ ;x )= (l-x) m 2 F 1 [ ; , 2i 

\ — m / V — w a; — 1 / 
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which is the limit case a := —m, b := n + 1, c — > — m of Pfaff's transformation formula 

/ a,b \ . , . „ / a, c — 6 x \ , s 

'HiT 11 " 1 ' iF, ( • ; ^t)' (2 ' 9) 

see |H (2.2.6)]. 

Remark 2.2. For the proof of (jl.ip from (|2.3[) as given above, one may start with a weaker 
form of (12. 3|) as follows: 

1 = (1 _ x) n+l q m , n {x) + x m+1 r m ,„(x), (2.10) 

where q m ,n{%) an d r m ^ n {x) are polynomials of degree < m respectively < n, so not yet necessarily 
explicitly given. Since (1— x) n+1 and x rn+l are polynomials without common zeros of degree n+1 
respectively m + 1, we can recognize (|2.10p as a Bezout identity, where q m ,n{x) and r m)Tl (x) will 
uniquely exist as polynomials of precise degree m respectively n (see for instance |5 5 Theorem 
6.1.1]). It was in this way that Daubechies §6.1], in the symmetric case m = n, proved (jl.ip . 
Also note that the symmetry of (I2,10p together with the uniqueness and degree properties of 
Qm,n(x) and r m>n (x) already imply that qm,n{x) = r n ^ m (l — x), without explicit computation. 

Equivalent to the Bezout identity approach, (|2.10p can be seen as a partial fraction decom- 
position 

1 Qm,n\X) fm,n\Xj 

x m+l(l _ x )n+l ~ x m+l (1 - X ) n+l 

with q m>n {x) and r mjn [x) of degree < m respectively < n, cf. Jagers' proof of (jl.ip in [16]. 

3 A proof by induction 

The following proof by induction was essentially given earlier by Kleiman [15] and, for m = n, 
by Daubechies [U Lemma 4.4]. We have to prove (|1.2p . with p m>n (x) given by (|1.3p . First note 
that (jl.2p holds for n = 0, and hence, by symmetry, also for m = 0. Indeed, 

m i _ _m+l 



Pm,o( x ) = (1 - ^) V X fc = (1 - X) — = 1 

1 — x 



k=0 

and po,m(l - x) = x m+1 , so p m ,o(a:) +Po,m(l - x) = 1. 

Now we prove (jl.2p by induction with respect to m + n. We already saw that it holds for 
m + n = 0, i.e., for m = n = 0. Now suppose that (jl.2p holds for all m, n with m + n = N — 1. 
Let m + n = N. Then we already proved (jl.2p if m = or n = 0, so we may assume that 
m, n > 0. Substitute the recurrence relation 

n + k\ fn+k — 1\ fn + k — 1 
J = VA;-1 J + V A; 

for binomial coefficients into (11.31). Then 



p m , n (x) = (i - xr 1 z( n t!i V +(!- -) n+1 E ( n + £ 

fc=l ^ ' fc=0 ^ 
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= x(l-x) n+1 J2 [i Jx l + (l-x) Pm , n ^(x) 

= Xp m -i tn (x) + (1 - X)p m)n -l(x). (3.1) 

Hence 

Pn,m(l ~ x) = xp njm -i(l - x) + (1 - x)p n -l,m(l ~ (3-2) 
Adding (|3TT|) and (pT2j) gives 

Pm,n(l) +Pn,m(l - x) = x(Pm-l,n(«) + Pn,m-l(l ~ %)) + (1 ~ ^)(Pm,n-l(^) +Pn-l,m(l ~ z)) 
= X + (1 — x) = 1 



by induction. This completes the proof of f) 1 . 2 j) . 



4 A proof by repeated differentiation 



Here we give a proof of JTT6J) which was sketched by Pieter de Jong in an earlier version of |10| , 
First note that the case n = of (11. 6p is evident. It is essentially the summation formula for 
the terminating geometric series: 

i T m+1 

— = l + x + x 2 + ... + x ™ + — . ( 41 ) 

1 — x 1 — X 

Now we can prove (|1.6p by induction with respect to n (for each n for general m). For n = 
we have (|4.ip . which is evident. Apply the operator (n + l)" 1 d/dx to both sides of (|1.6|) . The 
left-hand side becomes (1 — x) _ ™~ 2 , the right-hand side becomes 

( ra + 2 )fc Jfc , ■m V- m + ! + _ ^k-n-l , ™+l V n ~ + 1 ( m+1 )fc n _^fc-n-2 

2^ fcj X+X Z^ n+ 1 fc J ^ ^ +x 2^ n+ l k \ y y x > 

fc=0 k=0 k=0 

So we will have proved (|1.6p with n replaced by n + 1 and m replaced by m — 1 if we can show 
that 

x myrn + l (m + l) fc _ fe _ n _ x +1 ^ n-fe + 1 (m+l)» _ fc _ n _ 2 
^n + l k\ y 1 ^ n + l k\ y 

k=0 k=0 

is equal to 

n+l / % 

x mJ2^(l-x) k - n - 2 . (4.3) 

k=0 

In order to show this, rewrite x m+1 in the second term of (j4.2j) as x m+1 = x m — x m (l — x), by 
which (14.21) becomes 



m ^ m-n + fc (m + l)fc , ,k-n-i , m ^ n-k + l (m + l) fc . fc _ n _ 2 
2- n + l fc! 1 j + 2- n + l fc! 

fe=0 fe=0 



6 



^ m m-n + k-l (m + l) fc _i _ x)k -n-2 + x ™ y- n-k + l (m + l) fc . _ xfc_„- 2 

^ ra + 1 (fc-1)! ^ n + 1 fc! 

fe=l v ' fc=0 



= x m (i - x)- n ~ 2 + x m (! - *O fc ~ n ~ 2 + ^^x m (l - x)-\ 

k=l ' v n + J- 

which equals (14. 3p . This completes the induction step. 



5 A proof by generating functions 

In [17] a proof by generating functions was given for the n- variable generalization (|1.5p of (11. ip . 
Of course this proof can be specialized to a proof by generating functions of (jl.ip . Such a proof 
of (jl.ip was also communicated to us by Helmut Prodinger, independently from |17| . Because 
the one-variable case is more simple, we give the proof here. 
Fix x e (0, 1). For u, v e (0, 1) let 

/(«,«;*):= 2 ^ n (l-^ +1 E( n ^)^ (5-1) 

m,n>0 k=0 ^ ' 

Then 

f(v,u;l-x) = ^ u m v n x m+1 Y,(^ )a-x) k - (5-2) 

m,n>0 k=0 ^ ' 

From (|5.ip we have 

/(«, «;*) = -^-2> b (i- *) n+1 T,( n+ k k ) m* 

U n>0 fc>0 ^ ' 



1 - u ^ 



- 



(1 - ux) n+1 

n>0 v ; 

1 -X 1 



(1 - u)(l -ux) i _ ?!(lzg) 

1 -X 1 



Hence 



/(v,u; l-x) 



1 — u 1 — ux — v(l — x) 

X 1 



1 — v 1 — ux — v(l — x) 



and 



f(u,v;x) + f(v,u; 1 - x) = 3— ( - — - + 



1 — ux — v(l — x) Vl — u 1 — v J (1 — u)(l — v) 

So _^ 

f(u,v;x) +f(v,u;l - x) = ^ u m v n , 

m,n>0 

and combined with (|5.ip . (|5.2I) this yields (jl.ip by taking the coefficient of u m v n . 
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6 A proof by weighted lattice paths 



Consider all lattice paths from (0,0) to (m + l,n+ 1) in the planar integer lattice (using only 
unit east and north steps). Such a path P consists of m + n + 2 successive unit steps s k {P) 
(k = 1, 2, . . . , m + n + 2). Let be the path P terminated after k steps. The weight w(P) of 
a path P is defined to be the product of the weight of the respective steps s of the path, i.e., 
w(P) = Y\ se pw(s). Define the weight function w as follows: 



w{(i,j) -> (i + 1, j)) :-- 
-> + 1)) := 



x (j<n + l), 

x + y (j=n + l), 

y (i < m + 1), 

x + y (i = m + 1). 



Since for each k G {1, . . . , m + n + 2} and for each truncated path Pfc-i we have 

^2 w{s k {P k )) = x + y, 

Pf, Pk\s k (P k )=Pk-i 

we find by induction that Y^p k w {Pk) = {% + y) k , and hence 

Y J w(P) = (x + y) m+n+2 . (6.1) 
p 

On the other hand each path P ends either with a vertical step or with a horizontal step. 
Consider first the paths which end with a vertical step. Then the last horizontal step will 
be (m,k) — > (m + 1, fc) for some k G {0,1,... , n}. For given all such paths have weight 
x m+l y k {x + yy i ~ k + 1 and the number of such paths is ( m k k )- Hence the sum of the weights of 
all paths which end with a vertical step equals 

x m+i ib( m t k ) ykix+y)n ~ k+1 - (6 - 2) 

fc=0 ^ ' 

Similarly, the sum of the weights of all paths which end with a horizontal step equals 



m i 
h—n N 



n + k \ X \ X + y) m - k +\ (6.3) 



, k 

k=0 x 

Since (16. ip = (16. 2H + (I6.3p . we have obtained (|2.6p with both sides multiplied by x + y. 

Remark 6.1. For < x < 1 and y := 1 — x we can give a probabilistic interpretation of the 
above results. Now consider all lattice paths from (0, 0) to (m + 1, n + 1), where each following 
step has probability 1 if there is only one possible step, and otherwise probability x if the step is 
horizontal and 1 — x if it is vertical. Then (|6.2p gives the probability that the last step is vertical 
and (|6.3p the probability that the last step is horizontal, and the sum of both probabilities is 
necessarily 1. Thus we have a probabilistic proof of (|l.lj) . 
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Many probabilistic proofs of (jl.ip were earlier given, see [16], [18] . [22] . They are essentially 
all equivalent to the one given in the previous paragraph. Zeilberger's [22] proof (phrased by 
him for m = n) is particularly succinct. For general m, n it reads as follows: 

Toss a coin (with Pr(head) = x) until reaching m + 1 heads or n + 1 tails. Then equate the 
probability, 1, of finishing after at most m + n + 1 tossings with the sum of the probabilities of 
all the final outcomes. This yields (jl.lj) . 



7 A proof using the beta integral 

By the evaluation of the beta integral (see [U §1.1]) we have for x G (0, 1): 
(m + ra+1)! f 1 



ml n! 

n -+ 
m! n! 
(m + n 
m! 



/ t m (l-i) n cft 
Jo 

(m+n+1)! /%- ( i-trdt+ (m+ ;+ 1)! fni-trdt 

Jo m\ n\ J x 

n -±^ r^(i-r^+ (m+ "t 1)! r x m-tTdt. (7.i) 

!n! y m!n! J 

Then (jl.ip will follow from (|7.ip if we can prove that 

(m + "+ 1)! r ^ - *)" * = Y " (7-2) 

m!n! V ^ fc! 

But (|7.2p follows by the string of equalities 

r t m (l - t) n dt = x m+1 f 1 s m (1 - s + s(l - x)) n ds 
io io 

- ^ t (I) (i - [ »™+*u - .r* * - ^ " * 

fc=0 v 7 - 70 fe=0 

The integral on the left-hand side of (17. 2|) is an incomplete beta function, which is usually 
expressed as a hypergeometric function (jl.7p (see §2.5.3, p. 87], also for m, n complex with 
Rem > — 1): 

B^m + l,n + l):=[\ m (l-t) n dt = -}—x m+1 2 F 1 (~ n ' m + 1 -,x) (x G (0, 1)). (7.3) 
Jo m+1 \ m + 2 J 

The proof of (|T.3|) is by binomial expansion of (1 — t) n . Then (|7.ip takes the form 

T(n + m + 2) . , Tin + m + 2) 

1 = — t—, ^ BJm + 1, n + 1 + -— i — i— ■ fli_a,(n + 1. m + 1) 

r(m + l)r(n + 1) 1 ; r(m + l)r(n + l) v y 

(x G (0, 1), m,n G C, Rem, Ren > -1). (7.4) 
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The right-hand side of (|7.2[) can be written as a terminating hypergeometric series (j 1 . 8 f) . 
Then combination of (|7.2p and (|7.3p yields 



^ ( " n ' m + 1 ; l - x ) = - m+1 ^ ( 7 X 1 ^ ) • ( 7 - 5 ) 

\ — n J (m + l)!n! \ m + 2 J 

Alternatively, (|7.5p can be proved as the limit case for c — ► — n of Pfaff's identity 

f—n,b \ (c — b) n ( —n,b \ 

oF\ • 1 — x = — ^ a-fi , ; x ( n nonnegative integer), (7.6) 

V c / (c) n V&-c-n+l/ 

see [U (2.3.14)]. 



8 Extension of the identity to non-integer m and n 

By (J72D, CL3]) the formula 

TCttj -)- fi -\- 2) 

p m ,n(x) = r ( m + ^ r ^ n + 1 ) B i-x{n + 1, m + 1) (x G (0,1), m,n G C, Ren > -1) (8.1) 

extends (|1.3|) to non-integer values of m, n. Then, by (|7.4p . the identity (jl.2j) holds for x G (0, 1) 
and m,n G C with Rem, Ren > —1 if p m ,n(x) 1S given by (|8.ip . Moreover, by Carlson's theorem 
(see for instance Titchmarsh |20|. §5.81]) this is the unique extension 

Pm,n(x) + q m , n (x) = 1 (x G (0, 1), m, n G C, Rem, Ren > -1) 

of (|1.2p such that for nonnegative integer m,n we have p m .n(x) = <Zn,m(l — x) given by (|7.2p . 
Pm,n{x) and q m ,n{x) are analytic in m, n for Rem, Ren > —1 with x fixed, and p m ,n{x) and 
Qm,n(x) satisfy, for some c G (0,7r), estimates 0(e c ' m l) and 0(e c ' n l) as Rem, Ren > 0. Indeed, 
fix m with Rem > —1 and x G (0, 1). Then, in the right-hand side of (|8.ip we have for Ren > 0: 

1—x _ (1 — £)R em +l 



\B 1 - X (n + l,m + 1)| < 5i_ x .(Ren+ l,Rem + 1) < / (1 - t) Rem dt 

JO 



Rem + 1 

and (as a consequence of the asymptotic formula for T(z), see [TJ Theorem 1.4.1]) 



T(m + n + 2) 



r(m+ l)T(n + 1) 



0(\n 



Rem+h 



Hence, for x,m fixed as above, the right-hand side of (|8.ip is 0(e c ' n ') as Ren > for arbitrary 
small c > 0. We can estimate the other cases in a similar way. 
Alternatively, we may write (jl.3p as 
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Consider (|8.2p for x £ (0,1) and n E C with Ren > —1, and then try on it the fractional 
extension of finite sums proposed by Miiller & Schleicher [11], [12] • Since for £ C with 
Re fc > we have 

/(*) : = r( ^ ( " t) r(t + 1) xk{1 ~ xT+1 = o(1) asRe ^°°' 

their recipe of fractional extension (see [HI (10)], [HJ top of p. 5]) of the sum J2T=of( k ) * s 

oo 

) = Y,(f( k )-f( k + m + 1 )) 

k=0 

= (1 - ^ ( n ± ^ _ (1 _ x )n+l ^r(n + m + fc + 2) 



A _ n M v y ^r(n + l)r(m + A; + 2) 

T(n + m + 2) ^ _ +1 /n + m + 2,1 \ 

T(m + 2)r(n + l) 1 j 2 \ ^ + 2 ' J 

T(n + m + 2) m+1 ^ f-n,m + l 



= 1 

r(m + 2)r(n + l) V m + 2 

T(n + m + 2) 

= l~ r / , 1UV m+l,n + l . 8.3 

1 (m + 1)1 (n + 1) 

Here we used Euler's transformation formula [1, (2.2.7)] and (|7,3p . Note that for x € (0,1) 
and m, n G C with Rem, Ren > —1 the extension of p m ,n{x) defined by (|8.ip is equal to the 
extension defined by (|8.3j) . This equality is given by (|7.4p . Curiously, this equality is also the 
extension of the identity p. 21) . 

9 Three-term hyper geometric identities 

We can write (ll.6p as 

u 3 = u 1 + u 2 , (9.1) 

where 

„ lW:= £fe^ = 2fl (-^+V), (9 . 2) 



fc=0 

„.,,,., ._ ,-1 y + _ x )fc-n-l = x m+l (1 _ x) -n-l 2 F X f ™ + 1 ; 1 - x 

fro fc! V " n 



(m + 1)! n! \ m + 2 / 

u 3 (x) := (1 - x)-"- 1 = (1 - x)-"- 1 2 Fx f °' " m " n " 1 ; 1 - xV (9.4) 



Here the third identity in (|9,3|) is (|7.5|) . Now ui, ii2 and 113 are three different solutions of the 
hypergeometric differential equation 

x(l — x)u"{x) — ((n + 2)x + m(l — x))u'(x) + m(n + l)u(x) = 0. (9-5) 
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Indeed, consider Rummer's 24 solutions of the hypergeometric equation 



z(l - z)u"(z) + (c-{a + b+ l)z)u'(z) - abu(z) = (9.6) 

in §2.9] with (a,b,c) := (— m,n + 1, — m). Then Ui,U2, U3 above (uz up to a constant factor) 
are equal to (1), (17), (21), respectively, in §2.9]. However, this is for u\ and 113 only a 
formal proof, because there occurs a lower parameter in the hypergeometric function which is 
a nonpositive integer. For a rigorous argument for u\, consider the solution 2.9(1)] of (|9.6p 
first for (a,b,c) := (— m, n + l,c) and then let c — > — m. Also, for U3, consider the solution 
2.9(21)] of (|9.6p first for c, o := — b — m — n — 1 with 6 general, and then let b — » n + 1. 

For the general theory of solving the hypergeometric differential equation (|9.6p see [1, §2.3]. 
In general, for fixed a, b, c, and on a simply connected domain in C which avoids the singular 
points 0, 1 (and 00), one can choose two linearly independent solutions and have the general 
solution as an arbitrary linear combination of these two solutions. For a particular solution 
the coefficients in the linear combination can be found from (possibly asymptotic) values of 
the solution at two of the three singular points. In our case of solutions ui, 112, U3, given by 
(|9.2p - (|9.4p . the solutions are rational, so they exist as one- valued functions on C (possibly with 
a pole in 1). If we would a priori know only that us(x) = Au\[x) + Buzfa) then we can compute 
A = 1 by putting x = and we can compute B = 1 by multiplying both sides of the equality 
by (1 — x) n+1 , next putting x = 1, and then using the Chu-Vandermonde identity [1, Corollary 
2.2.3] for the evaluation of 2F\(m + 1, — n; m + 2; 1). 

The case discussed here is the case m = in Vidunas [211 §8] (trivial monodromy group). 
In this way he obtained (II. 6p as the case m = of the identity at the end of his section 8. 

For a > and n, m nonnegative integers we will now prove the following more general 
three-term identity: 

-n-l/1 ~ (m + l,-a 



(1-^-^(1-^)-%^ ' ' ;* 

n + m + 2 



(n + lWl z m+l( 1 _ z yn-l( 1 _ z -lya 2 pf n,m+l^ 



(n + a + l) m+ i V —n — at 

(m + l)„+i „ / —m, n + 1 \ / „. . , » 
(m + a + l)n+l \ -m — a J v 

This formula makes good sense on the indicated domain, since 2-^1 (&> b; c; z), originally defined 
as a power series for \z\ < 1, has a unique analytic continuation to C\[l,oo). Hence, the 2F1 
on the left is uniquely defined for z £ [0, 1]. Also, (1 — z~ l ) a is uniquely defined for z £ [0, 1]. 
The two 2-P1 's on the right, being polynomials, are defined for all z G C. If one wishes, one may 
rewrite (1 - z)"™"^! - z" 1 )" as (-l) n+1 z~ n ~ l (1 - z" 1 )- -™- 1 . 



For a —t the identity (|9.7p tends to the identity (|1.6p . In fact, (|9.7p . which is of the form 
U3 = u\ + U2 (see (19.1j) ). has the terms u\, 112 and 113 as solutions of the differential equation 

z(l - z)u"(z) - ((n + 2)z + m(l - «) + a)-u'(z) + m(n + l)u(z) = 0, (9.8) 

i.e., the hypergeometric differential equation (|9.6p for (a,b,c) = (—m,n + 1,—m — a). Indeed, 
see the solutions 2.9 (18),(1),(14)] of ([23]). This gives «i (after substituting (06])), n 2 and 
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us, respectively. For a | 0, (|9.8p tends to (|9.5p . and the solutions ui,U2,us of (|9.8p tend to the 
solutions ui,U2 5 ii3 of (|9.5p . 

The case of (|9.7p that a is a nonnegative integer is essentially the general case of the identity 
at the end of section 8 in Vidtlnas [21] . Just transform the 2-F1 on the left of (|9.7p by first reversing 
the order of summation in the hyper geometric series and next applying Pfaff 's transformation 
formula PH . 

For the proof of (|9.7p , start with the three-term identity 

/m + 1, — a \ _ T(n + m + 2)r(n + a + 1) _ m _! /m + 1, — n _ 1 \ 
2 1 Vn + m + 2 ' V ~ T(n + l)r(n + m + a + 2) Z 21 ^ -n - a ' ~ 2 7 

, r(n + m + 2)r(-n - a - 1) ^n+i+a v ( n + ^- m ., v -i 

T(m + l)r(— a) \ n + a + 2 

(z£{0}U[l } oo)), (9.9) 

see 2.10(4)] (or [H (2.3.11)] combined with (|2~3|0 . Note that we don't have to exclude 
z G (—oo,0) in (19. 9p because the two 2-^1 's on the right are terminating. By (|7.6f) the last 
hypergeometric function on the right can be replaced by 

(a + l) m (n + 1,-m _ t 

2*1 I \Z 



(n + a + 2) m \ — m — a 

In the identity which thus results from (|9.9|) . first replace 2; by z -1 and next multiply both sides 
by (1 - z)- n ~ l {l - z~ l )- a . This yields d577|l . 

10 A multivariable generalization 

A multivariable generalization (|1.5p of (jl.ip was proved in [17] by generating functions (see the 
one- variable case of this proof in section [5|) , while a probabilistic proof, immediately generalizing 
the one- variable case discussed in Remark 16. 1[ was indicated in [17J and [22] . We will give in 
section [T2] a different proof of (|1.5|) , which will generalize the proof of (|1.1|) in section [7] using 
the beta integral. 

Let us reformulate (|1.5j) in other notation, and let us also give this identity in homogeneous 
form. Let s := x\ + ■ ■ ■ + x n . Define 

ai a n -i 1 1 -1 \ 

f ( T , r,. \ ._ r an+l V V 1 ' " + 

Jai,...,a n ^l, • • • .- X n 2^, ■ ■ fci!...fc_ 1 ! 

fci=0 fc n _i=0 

x x^ 1 . . . x k n n _'i s a i+-+ a »-i-( fc i+-+ fc ™-0. (10.1) 

Then 

(XI + • • • + X n ) ai + - +a " +1 = ^ /a CT(1) ,...,a CT(n) (^(l), • • • , ^(n)), (10.2) 

a 

where summation is over all cyclic permutations a of 1, 2, . . . , n. For x\ + ■ ■ ■ + x n = 1 identity 
(fT02l simplifies to 

1 = 5^/o <r( i ) ,...,aa(„)( a; ff(l)>---> a; «r(»)) ( 10 - 3 ) 
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with 

a\ a n -i 



fai,...,a n (x\, ■ ■ ■ , X n ) — X n n + • • • ^_ , • • • x n -l ■ (10-4) 

fci=o fc„_i=o "' '"~ L " 



Identity (|10,3p is a reformulation of (jl.5p and (|10,2p is the homogeneous form of (|10.3p . 
Remark 10.1. For n = 2 we can write ()10.3p . (|10.4p as 



1 = /m,n(x, 1 " X) + f n , m (l ~ X, x) with f m , n ( X , 1 - x) = (1 - x) n+l ~ ^ ' ' 



fc=0 



So we have (jl.ip . The case n = 3 of (|10.3p , (|10.4[) is also noteworthy as a two- variable analogue 
of ([LI]): 

1 = fa,b,c( x ' VA ~ x ~ y) + fb,c,a(y, l-x-y,x) + / c , a ,&(l - x-y,x, y), (10.5) 
/ 0)6iC (x, y, 1 - x - y) = (1 - x - y) c+1 E (C+ kin +l ( ia6 ) 

k=0 1=0 

Identity (|10.5p is also given in [T7j. 

If we divide both sides of f)10.5|) by (1 — x — y) c+1 then the resulting identity has the form 

(1 - x - y)-^ 1 = f(x, y) + x a+1 g(x, y) + y b+1 h(x, y), 

where f(x,y) is the polynomial consisting of all terms of the power series of (1 — x — y) _c_1 
which have degree < a in x and degree < b in y, while g(x, y) and h(x, y) are power series in x 
and y. There does not seem to be an a priori symmetry argument which settles (I10.5P from this 
observation. For instance, if we try to imitate the proof for n = 2 in section [2] then we have to 
consider the homogeneous version of (|10.5p given by (110.2H for n = 3: 

(X + y + Z ) a+b+C+1 = f a ,b,c( X > y> Z ) + fb,c,a{y, Z, x) + fc,a,b{z, X, y). 

The problem is where to put on the right the terms "1k,l,mX k y l z m (k + l + m = a + b + c+l) in 
the expansion of (x + y + ^) a + b + c + 1 . Certainly, we can uniquely put all terms , yk,i,m xk y l z m with 
k < a, I < b in f aj b jC (x,y,z), all terms with I < b, m < c in fb,c.a(y, z, x), and all terms with 
m < c, k < a in f c ,a,b( z ,x,y), but there is no clear rule where to put a term in which only one 
of the three inequalities k < a, I <b, m < c holds. 

Remark 10.2. Note that we can formally express the right-hand side of (|10.6p in terms of 
Appell's hyper geometric function Fi (see [3 §5.7.1]): 

fa,b,c( x i y, 1 ~ x - y) = (I - x - y) c+1 F 2 (c + 1, -a, -b, -a, -b, x, y). 

However, due to the nonpositive integer bottom parameters we cannot transform this F 2 function 
by [3 5.11(8)] similarly as we transformed the 2 -Pi function in Remark 12. 1[ 
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Similarly to the case n = 3, the multisum on the right-hand side of (|10.4p can be formally 
written as a Lauricella hypergeometric function Fa (see [21 Ch. VII] or |13[ (8.6.1)]): 

fai ,...,a n (%X j • • • j Xn) = ^n™ A (fln 1 ; ^1 j • • • j Q"n— 1 j j • • • ; ^n— 1 > 3^1 > • • • j 2-n— l) 

(xi + --- + x n = l). (10.7) 

Here the .F^ 1 ^ has to be interpreted as 

lim F^ n ~ 1 ^(o n + l,-ai,...,-o n _i 3 6i,...,6n-i;aJi,--- J »n-i)- (10.8) 

(6i,...,6 n _i)-+(— ai,...,— a n _i) 

11 A PDE associated with the multivariable identity 

From the expression in (|10.7p of f ai ,...,a„( x i, ■ ■ ■ > ^n-i> 1 — sci — • • • — a? n _i) in terms of a Lauricella 
hypergeometric function, we will derive a PDE for f ai ,...,a n - Consider first the system of PDE's for 
thfii^ n-1) in (mm as given in [21 Ch. VII, §XXXIX] (for n = 3 we have Appell's hypergeometric 
function F<i and then the system of PDE's is also given in [7J 5.9(10)]). After taking the limit 
for . . . , b n -\) — > (— a\, . . . , — a n -\) we obtain that 

u(xi, . . . ,x n -i) := F A n l '(on + 1, -ax,. . . , -a n -i, -ax, • • • , -a n -i; a;i, - - • , £n-i) (H-l) 
satisfies the system of PDE's 

Xi(l - Xi)dju - Xi XjdjdiU - (en + (a n + 2)xi)diU + ^jdj ^ + (a n + l)ojti = 

(i = 1, . . . ,n- 1). 

Here di denotes d/dxi. The sum of the n — 1 PDE's equals 

y] Xi(l - Xi)dfu — ^ XiXjdidjU + (ax H h a n _i — a n — 2) Xj<9jU - a^u 

i i^j i i 

+ (a n + l)(ax + ■ ■ ■ + a n -x)u = 0. (11.2) 

Proposition 11.1. The function (jll.ip is £/ie unique solution, up to a constant factor, of (jll.2p 
which has the form 



ai ln-1 



n(xi,...,X n _i) = ^ ... ^ Ifci,...,*;™-!^ 1 •••^n-1 1 - ( 1L3 ) 
fci=0 fc n _i=0 

Proof Computation of the left-hand side of (jll.2p with u := x^ 1 . . . x^Zx yields 



{(ax - k X ) + • • • + K_i - fe,»-l)) K + 1 + fa + ■ • • + fen-l)xf • • • XnSi 

+ (E - ^ - *K rl J x l ■ ■ ■ x n-l- 
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It follows that u of the form ([TO]) satisfies (fTO|) iff 

((ai — ki) H h (a n _i - A;„-i)) (a n + I + k\ -\ h fc n -i)7*i,...,fe„-i 

n-i 

+ + ^ " a iHl,.~to+h~,kn-l = - ( 1L4 ) 

1=1 

Give some value to 7 ail ... o w _i- Then we see from ([11. 4p by downward induction with respect to 
ki + • • • + k n -i that all coefficients 7fc lv .. i fc„_ 1 with < k{ < (i = 1, . . . , n — 1) are uniquely 
determined by this initial value. 

In passing we see that (|11.4|) is satisfied by 

(a n + l)fei+-+fc„_i 



^>-,kn-i ■-- ... k n -\\ ' 

Thus we have also proved from scratch that u given by (|11.1|) satisfies (|11.2p . □ 
By some computation, we see that 

v(xi, x n -i) := (1 - xi x n _i) an+l u(xi, x„_i) 

— fai,...,a n ("^l) ■ ■ ■ j %n— 1; 1 -^1 ' ' ' ^n- 1) 

satisfies the PDE 

n— 1 n— 1 

y^acj(l - Xi)dfv - 2~^2xiXjdidjV + yY(ai H h - aAdiV = 0. (H-5) 

i=l i<j i=l 

Clearly, the function v := 1 satisfies ([11. 5|) . Furthermore, by straightforward computations we 
see: If v satisfies .5[) then the function (xi, . . . ,x n -\) t— » v(l — x\ — ■ ■ ■ — x n -i,X2, ■ ■ ■ ,x n -i) 
satisfies (|1 1.5[) with a\ and a n interchanged. Thus we have proved: 

Theorem 11.2. For all permutations a of 1,2, ... ,n the functions 

(xi, . . . ,X n -i) l ^ /o ff(1) ,...,o . (n) (a;£r(l))-- -,x a (n)) (xi H hx n = 1) 

are solutions of ([11.5[) . t/p to a constant factor they are the unique solutions of (jll.5p o/ the 
form 

a cr(l) a a(n-l) 

v(x u ...,x n . 1 ) = xl^ ) )+1 ■■■ Yl Tfe!,..,^-!^!)---^-!) (x 1 +---+x„ = l). □ 

fci=0 fc n _i=0 



Now consider a solution v(xi, . . . , x n -i) of (111.51) . let x n be a variable independent of 
x%, . . . , x n -\, and let <f> be an arbitrary function of that new variable. Then trivially (jll.5p holds 
with v replaced by 4>(x n )v(xi, . . . , x n -\). Now pass in this PDE to new variables y\, . . . , y n by 

_ m _ Vn-1 _ 

x\ — , • • • , x n — i , x n yi + ■ • • t y n , 

yi H h y n yx H 1- ?/n 
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or equivalently, 

Hi — ^iX n , . . . , Un—l — X n —iX n , y n — (1 X\ ■ x n _i)x n . 

Then we obtain that the function w(yi, ...,y n ) ■= 4>{yi H + Vn) v( yi+ n , +yn , ^+",7+^ ) 

satisfies the PDE 

n n 

^yi(yi-\ VVn- Vi)dfw -2y~] yiyjdidjw + ^ Nai H h a n )yi - a»(yi H h y n )) <9i™ = 0, 

i=l i<j i=l 

(11.6) 

where <9j denotes d/dyi. Thus by (jlO.ip . (|10.4p and Theorem 1 1 1 . 2 1 we have proved in particular: 

Theorem 11.3. The function io(yi, . . . ,y n ) := / ai ,...,a„(yi> ■ ■ ■ ,Un), defined by (jlU.ip . satisfies 
P1-6P. Similarly, the functions w(yi, ... ,y„) := /o ct(1) ,..., Mre) (^(l), • • • , yo-(n)) sa^s/y flll-6|), 
where a is a permutation o/l, . . . , n. ylfeo, the function w(yi, . . . , y„) := (yi + • • •+y n ) ai ~ l han+1 
satisfies (jll.6p . 



12 Splitting up Dirichlet's multivariable beta integral 



Just as (|1.6p can be obtained by splitting a beta integral into two parts and evaluating the 
resulting incomplete beta integrals, we can prove and interprete f)10.3|) by splitting Dirichlet's 
(n — l)-dimensional beta integral with nonnegative integer exponents into n parts. For conve- 
nience, we will work here with an n-dimensional beta integral. 

Let A n be the simplex in R n wich has vertices and the standard basis vectors e\, . . . , e n . 
Let ai, . . . , a n , b be complex numbers with real part > —1. Then Dirichlet's integral is as follows. 

I ai ,..,a n+1 := [ t? • ■ ■ C (1 - h tnT^ dt x ... dt n = ^ Ql + 1 ^- r(a ; +1 + \ ) , (12.1) 

Ja„ r(ai h a n+ i + n + 1) 

see [TJ Theorem 1.8.6] or [6i, Exercise 7.2.6] for a straightforward proof, and [8] for its history. 
Note that I ai ,...,a n +i * s symmetric in a±, . . . 

Now take x = (x\, . . . , x n ) within A ra and let An\x) (i = 1, . . . ,n + 1) denote the simplex 
in lR n which has a vertex x and n further vertices taken from 0, ej, . . . , e n where ei is deleted if 
% = 1, . . . , n and is deleted if i = n + 1. Require that ai, . . . , a n+ i are nonnegative integers. 
Define 

41., an+ »:= / n ^...C(l-*i tn^rfti.-.c^. (12.2) 

For any (yi, . . . ,y n ) € W l put y n +i := 1 — yi — • • • — y n - Then, for any permutation a of 
1,2, ...,n+ 1 (i.e., cr € SWfl), the map (yi,...,y n ) (y<x(i), • • • , y<r(n)) is a diffeomorphism of 
A n with Jacobian having absolute value 1. Thus we have the identity 

Ial...,a n+1 ( X ^ ■ ■ ■ > x n) = 4^, , W ?a CT(n+1) (^(1) , • • • , X*(n)) G S„+l). (12.3) 

Define 

fa 1M (*) := (1 - *i x n f^ £ £ fcl l + \' +kn ^ • • • ^ ' (12 ' 4) 

fc 1= fc„=0 1- • • • n- 
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i.e., (jlO.ip with n replaced by n + 1 and with x n +i := 1 — x\ x n omitted in the argument. 

We have the symmetry 

fa 1 ,...,a n+1 \ x li . . . ,X n ) = fa tT(1) ,...,a tT ( n ),a n+1 {x cT (l)-, ■ ■ ■ , x a(n)) { a G ^n)- (12.5) 

Proposition 12.1. For nonnegative integers ai,... ,a n+ \ and for x within A n we have: 



I (n) (t) 



fa n +i,ai,...,a n (1 x l x ni x li ■ ■ ■ , X n — i). (12.6) 



L ai,...,a n +i 

Proof For convenience put b := a n+ \ and x' := (x±, . . . ,x n -\). Then 
I (n) (x) - 

C ( / t? 1 • • • C-l 1 (1-*1 tut dt X ... dt n _! J dt n 

+(1— (* n /a; n ))A„_i / 

x a n+ l f s a n f f t a, ^-i ^ _ ^ ^ _ ^6 ^ ^ j ds 

JO \Jsx' +(l-s)A n - 1 J 

x a n+ l j s a n (1 _ s) n-l I j ^ + (1 _ ^^cw _ _ _ + (1 _ s)^)-*-! 

J \M„_i 

.6 



x (s(l - xi x n ) + (1 - s)(l - si s n _i)) dsi . . . ds n _ij (is 



^1/ \kn—lj \lj 
x ^ s a„+feiH hfcn-l+/ ^ _ s ^aiH \-a n -i+b+n- i-fci k n -i-l ^ 

ST""' ■ ■ ■ C-r^ 1 (1 - S l S n-l) 6 "' ds l • • • 



fci=0 fc„-i=0 Z=0 

»1 

„a„+fciH hfcn-l+/ / 1 _ „\aiH han-i+t>+n— 1— fci k n -i-l 







'A n _i 

e EGO-drOCK-ei'd-.--^ 

T(a n + fci H h A; n _i + Z + l)r(ai H h a n _i + 6 + n — &i /c n _i - Z) 

T(ai H \- a n + b + n+l) 

r(ai - fei + 1) . . . r(a TO _ 1 - fc n _x + i)r(6 - z + i) 

r(oi H h a n _i + Z> + n — A;i A; n _i - Z) 

ai! . . .a n \b\ 



[fix H 1- a n + 6 + n)! 

ai an-l b 



1 V V V (°* + 1 )fel+- + fen-l+i fei kn-lfl ... Nl 

fci=0 fc„_i=0 i=0 

Iai,...,a n ,b fb,ai,...,an-i,a n 0- ~ X ~L ~ ' ' ' ~ x ni x li ■ ■ ■ i x n-l)- d 
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Theorem 12.2. Let a±, . . . ,a n +i an d x as before. Let i G {1,2, . . . ,n + 1} and let a be the 

cyclic permutation of 1, . . . ,n + 1 which sends n to i. Then 

- t ai,...,a n +i V^y / \ /-.n y\ 

r ~~ /a CT (n+l),a CT (i).---,a CT (n)V X CT(n+l)> ^(7(1)1 • • • ,X cr ( n _ 1 )J. \ lz -<) 

J ai,...,a n +i 

Proof By we have J °j ) "-°"+ 1 (x) = J Md--M"+d ! ! ! ! Now apply (j^. 

□ 

We have the obvious identity 

1 = X) °p S " + (x G A n , a!,...,a n+ i GC, Reai,...,Rea n+ i > -1). (12.8) 

By Theorem 112.21 this is for nonnegative integers a\,... ,a n +i equivalent with (I10.3h (with n 
replaced by n + 1 and with the functions f ai ,...,a n+1 defined by (|1U.4|) ). In the general case of 
([12. 8p we get an extension of f)10.3j) for non-integer ai, ... ,a n +i, just as we discussed in the 
one- variable case in section [5J The uniqueness of the extension if the terms satisfy estimates as 
in Carlson's theorem, as discussed there, also holds here. 

Remark 12.3. It is an interesting question (but for us a nontrivial open problem) to find an 
elegant looking evaluation of (|12.2j) which is valid for all complex a\, . . . , a ra +i with real part 
> —1, and which would generalize the evaluation (|T.3[) of the incomplete beta function. This 
would also give an n- variable generalization of (|7.5j) . i.e., of a limit case of Pfaff's identity (17. 6ft . 
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